THIS  REPORT  HAS  BEEN  DECLASSIFIED 
AND  CLEARED  FOR  PUBLIC  RELEASE. 


DISTRIBUTION  A 
APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


UNCLASSIFIED 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 


CAMERON  STATION  UEXANDRIA.  VIRGINIA 


RATED ja  3 YEAR  INTERVALS: 
DECLASSIFIED  AFTER  12  YEARS 
DOD  DIR  5200  10 


UNCLASSIFIED 


1 


Ql 


QW 


STRESS  CONCENTRATION  .DUE  TO  A HEMISPHERICAL  PIT 
AT  A FREE  SURFACE 


B>  A.  Eubanks 


A Technical  Report  to  the 

Office  of  Naval  Research 
Department  of  the  Navy 
Washington,  D.  C. 


Contract  N7onr-32906 
Project  No,  NR  035-302 


Department  of  Mechanics 
Illinois  Institute  of  Technology 
Chicago,  Illinois 

January  1,  1953 


STRESS  CONCENTRATION  DUE  TO  A HEMISPHERICAL  PIT 
AT  A FREE  SURFACE 

By 

R.  A.  Eubanks 


A Technical  Report  to  the 


Office  of  Naval  Research 
juepartment  of  the  Navy 
Washington,  D.  C. 


Contract  N7onr-32906 
Project  No.  NR  035-302 


Department  of  Mechanics 
Illinois  Institute  of  Technology 
Chicago,  Illinois 

January  1,  1953 


hXl  - n<\oL  . 


SUMMARY 


This  paper  contains  a solution  in  series  form  for  the  stresses  and 
displacements  around  a hemispherical  pit  at  a free  surface  of  an  elastic 
body.  The  problem  is  idealized  by  considering  a semi-infinite  medium 
which  is  otherwise  bounded  by  a plane.  At  infinity  the  body  is  assumed 
to  be  in  a state  of  plane  hydrostatic  tension  perpendicular  to  the  axis 
of  synmetry  of  the  pit.  The  present  method  of  solution  may  be  general- 
ized to  loadings  which  are  not  rotationally  symmetric.  Numerical  results 
are  given  for  the  variation  along  the  axis  of  symmetry  of  the  normal 
stress  which  is  parallel  to  the  tractions  at  infinity;  these  results  are 
compared  with  the  known  corresponding  numerical  values  appropriate  to  the 
two-dimensional  analog  of  the  present  problem. 


TABLE  OF  CONTENTS 


Page 

iv 


PREFACE  . « « 

LIST  OF  TABUSS ▼ 

LIST  OF  ILLUSTRATIONS ▼ 

SECTION 

I.  INTRODUCTION,  STATEMENT  OF  PROBLEM 1 

II.  PARTICULAR  SOLUTIONS  IN  SPHERICAL  COORDINATES 5 

III.  SOLUTION  OF  PROBIEM 13 

IV.  STRESS  EVALUATIONS,  ACCELERATION  OF  CONVERGENCE, 

NUMERICAL  RESULTS 17 

BIBLIOGRAPHY 2h 


i 


I 


iii 


J. 


PREFACE 

The  author  is  greatly  indebted  to  Dr.  E.  Sternberg  for  invaluable 
advice  and  guidance  in  this  investigation.  Thanks  are  also  due  to 
Dr.  M.  A.  Sadowsky  for  hrlpful  suggestions.  The  assistance  of  Messrs. 
S.  Parter,  K.  Fukuda,  and  W.  F.  Darsow  in  voluminous  numerical  compu- 
tations is  gratefully  acknowledged. 


iv 


LIST  OF  TABIES 


Table  Page 

1*  Numerical  Values  of  for  >>  = l/L 17 

2.  The  Stress  for  p = 1,  V = l/U 22 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1.  Meridional  Cross  Section  of  Pitted  Region, 

Spherical  Coordinates  • 2 

2 • Decay  of  the  Maximum  Stress  along  the  Axis 

of  Symmetry,  = <r^  = 1 at  Infinity. 23 


STRESS  CONCENTRATION  DUE  TO  A HEMISPHERICAL  PIT 
AT  A FREE  SURFACE 

I.  INTRODUCTION,  STATEMENT  OF  PROBLEM 


In  what  follows  we  investigate  the  stress  concentration  produced  by 
a hemispherical  pit  at  a free  surface  of  an  elastic  body.  The  problem  is 
idealized  by  considering  the  medium  to  occupy  a serai-infinite  region  con- 
taining a hemispherical  indentation  and  otherwise  bounded  by  a plane. 

With  reference  to  cartesian  coordinates  (x,  y,  z),  the  region  under  con- 

o o 2 V 2 

sideration  (Figure  1)  is  defined  by  z - 0,  r = (x4"  + y + z ) -1  if, 

for  the  sake  of  convenience,  we  assume  the  pit  to  have  a radius  of  unity. 

The  body  forces  are  assumed  to  vanish.  Sc  far  as  the  loading  is  con- 
cerned, we  confine  ourselves  to  the  axisymmetric  case  in  which  the 
tractions  at  infinity  constitute  a plane  hydrostatic  field  of  stress 
parallel  to  the  plane  z = 0.  It  should  be  emphasized,  however,  that 
the  method  of  solution  adopted  subsequently  may  be  extended  to  the  case 
in  which  the  plans  stress  field  at  infinity  is  no  longer  rotationally 
symmetric.  In  the  present  instance,  the  boundary  conditions  at  infinity 
assume  the  form, 


cr  , cr 
x*  y 


If  cr  , 7*  , Y , Y -*0  as  r 

’ z*  'xy'  'yz*  'zx 


00. 


(1) 


where  (o^,  ...  7^.,  ...)  aie  the  cartesian  components  of  stress.  With 
reference  to  spherical  coordinates  (r,  ©,  y),  introduced  by  the  mapping, 


x = r sin  © cos  y, 
y = r sin  © sin  y, 
z = r cos  0, 


0 - r < oo. 


o - s ' 7 r. 


0 - V - ? It- 

i ' 


(2) 


1 


3 


the  conditions  to  be  met  at  the  stress-free  plane  boundary  appear  as 


% c rre  * Yey=0 


(3) 


for  1 - r < oo,  & = tf/2,  0 - y - 2 7f , 


•whereas  a traction-free  surface  of  the  pit  requires 


<r 

r 


(U) 


for  r = 1,  0 - & * *72,  0^  2^, 

- 

in  which  (a,  ...,  7l_  , ...)  are  the  spherical  components  of  stress. 

The  two-dimensional  analog  of  this  problem,  that  is,  the  problem  pre- 
sented by  the  stress  concentration  around  a semi-circular  notch  in  a semi- 
infinite plate  which,  at  infinity,  is  in  a state  of  uniaxial  tension  paral- 
lel to  its  straight  edge,  has  received  repeated  attention,  kaunsell  £l^ 
obtained  a series  solution  to  this  plane  problem  by  first  extending  the 
uniform  stress  field  at  infinity  throughout  the  half-plane  and  then  remov- 
ing the  tractions  so  arising  on  the  notch  with  the  aid  of  a doubly  in- 
finite sequence  of  particular  solutions  of  the  field  equations.  This 
sequence  is  constructed  in  such  a manner  that  each  of  its  members  is 
singular  at  the  center  of  the  notch,  and  leaves  the  otraight  boundary,  as 
well  as  infinity,  free  from  tractions.  Subsequently,  Weinel  [^2]]  and 
C.  B.  Ling  [3],  using  bipolar  coordinates,  arrived  at  more  elegant  inte- 
gral representations  of  the  solution  to  the  plane  problem.  The  applica- 
tion to  the  present  space  problem  of  analogous  integral  approaches,  based 


1 


Numbers  in  brackets  refer  to  the  bibliograohy  at  the  end  of  this 


paper 


1 


on  the  use  of  toroidal  coordinates,  meets  with  great  analytical  difficul- 
2 

ties.  The  method  of  solution  adopted  here  may  be  regarded  as  the  three- 
dimensional  counterpart  of  that  employed  by  Uaunsell  £l^. 


One  is  led  to  a system  of  dual  integral  equations  the  kernels  of 
which  involve  the  Legendre  functions  of  the  first  kind  of  complex  index, 
P_l/2  + it ^)»  w*lere  the  index  parameter  t is  the  variable  of  integration. 
These  functions  have  not  been  investigated  extensively;  see  PH*  p.  li$l. 


I 

_J. 


II.  PARTICULAR  SOLUTIONS  IN  SPHERICAL  COORDINATES 


The  problem  characterized  by  the  boundary  conditions  (1),  (3),  and 

(U),  nay  be  attacked  by  any  one  of  the  various  available  stress-function 

3 

approaches  to  axisymmetric  problems  in  elasticity  theory.  le  shall 
utilize  the  approach  originated  by  Boussinesq  ^7^]>  according  to  nhich  the 
general  solution  of  the  displacement  equations  of  equilibrium,  in  case 
of  torsion-free  rotational  symmetry,  and  in  the  absence  of  body  forces, 
is  representable  as  the  sum  of  the  two  displacement  fields, 

20  [ux,  Uy,  uj  = grad  0,  (5) 

2G  [ux,  uy,  uj  = grad  (z  (£/)  - U(1  - p)  £o,  0,  (jj\,  (6) 


where 


V^Cr,  z)  *=  0, 


72  (j Hr , z)  = 0. 


(7) 


Here  u , u . u are  cartesian  conponents  of  displacement, 

^ ' * i h 

2 2 2 ^ ^ 

r = (x  + y' + z ) , Vc  is  the  Inplace  operator,  and  G and  p are 

the  shear  modulus  and  Poisson’s  ratio,  respectively. 

For  our  present  purpose*  it  is  expedient  to  refer  to  the  spherical 
coordinates  (r,  e,  y)  defined  by  Equations  (2).  The  basic  displacement 
fields  given  in  Equations  (5),  (6),  and  their  rssociated  fields  of  stress, 
were  transformed  into  general  orthogonal  axisymmetric  coordinates  in  £6]]. 
The  specific  forms  assumed  by  the  Boussinesq  solutions  in  spherical  co- 
ordinates appear  explicitly  in  pQ,  which  also  contains  the  particular 
solutions  obtained  upon  introduction  of  the  interior  and  exterior  spherical 
harmonics  of  integral  order  as  generating  stress  functions  (ft  and 

^ 

"See  [_5J,  Chapter  VIII,  and  for  references. 
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In  view  of  the  regularity  requirements  inherent  in  the  current  problem, 
we  limit  our  attention  to  the  two  aggregates  of  harmonic  stress  functions. 


0„(r,  0)  = r_n  P (cos  0) 


n 


(r,  G)  ~ r”n  P ( cos  0) 


n' 


(n  =■  0,  1,  2,  • • . ) , 


(8) 


in  which  PQ  denotes  the  Legendre  polynomial  of  order  n.  The  stress 
functions  defined  in  Equations  (8)  represent  exterior  spherical  harmonics 
and  give  rise  to  displacement  and  stress  fields  which  are  singular  at  the 
origin  r = 0,  are  otherwise  regular,  and  vanish  at  infinity.  The  par- 
ticular solutions  of  the  field  equations  generated  by  and  (p q will 

be  designated  by  C*J  and  [Cj,  respectively.  As  explained  in  Q0, 
considerable  simplifications  arise  from  a replacement  of  solution  QjJ 
with  the  linear  combination. 


OJ  = <2n  + 1)  [Cj  - (n  + 1*  - U t»  J. 

We  now  cite'’  the  displacement  and  stress  fields  appropriate  to 
solutions  pg  and  £bJ.  the  auxiliary  notations. 


(9Y 


p = cos  0, 

we  obtain  for  solution  OJ> 

(n  * 1)  P, 


p c sin  0, 


(10) 


20^  «=  - 


n 


n+  2 


2G. 


P Pn 
"n+7» 


(11)' 


4 In  Equation  (9),  aa  in  subsequent  equations,  the  letters  in  brackets 
denote  either  the  displacement  vector  field  or  the  stress  tensor  field,  and 
equality,  addition,  as  well  as  multiplication  by  a scalar,  are  to  be  inter- 
preted accordingly. 


p J ^ ^ 

LTOO  U^UaUAUIlO  J , \1|  ),  Ol  I OJ« 


dp_ 


The  argument  of  Pn  is  henceforth  understood  to  be  p and  P^  e 


1 


9 = (n,+  1)(y.,2i  P , 

r n+3  n* 

r 

% = 7^3  [Pn+  1 - <”  * + 2>  1>n]< 


VT^’ 

On  the  other  hand,  for  there  results 


n+  1 


_ n + 2 - 
rr®~  TvTJ  P Pn 


► 02) 


20u  «•  - (°+  p 

r n+  1 n+  1 

r 


(13) 


>£(n 


<r,  = (n  + 1)|  (n  + l)(n  + U)  - 2 >> 


1 Pn+  1 


% = ^?[(n+  D(“2-^  1 - 2 >>)  PQ+  i - (n-3  + UP)  P^] 
<ry=“FT?[(1"2>>)(n+1)(2n+1)  Pn+  i+(n*3+U>>)  P^J 


y m 


f^s  (n2  + 2n  - 1 + 2 pf)  p ?n+1 


n + 7* 


In  the  preceding  equations  u , u_  and  o* , 9"  , <r  , denote 

r o r q y ry 

spherical  components  of  displacement  and  stress,  respectively!  the  dis- 
placement of  My  and  the  stresses  7^,  7^ , which  vanish  identically 
by  virtue  of  rotational  symmetry,  have  been  omitted. 

The  physical  significance  of  the  sequences  of  solutions  C‘j>  og. 
(n  = 0,  1,  2,  ...),  was  examined  in  Reference  QQ.  Solution  C*0!]  m&y 
be  identified  as  that  appropriate  to  a center  of  compression  at  the 
origin,  whereas  is  the  solution  corresponding  to  a concentrated 

force  acting  at  r = 0 in  an  infinite  body,  in  the  direction  of  the 


< 


J. 


8 


7 

negative  t-axis. ' Moreover,  in  cartesian  coordinates, 

5?W- 

CbJ  ■ l^.ngn-  [Bj  ♦ <n  - 3 + U>>)  C*n- J> 

(n  = 1,  2,  3,  ...). 


(15) 


Hence,  these  two  aggregates  represent  sequences  of  successively  higher 
singularities  which  may  be  obtained  by  suitable  limit  processes  applied 
to  the  basic  axisymmetric  singularities  corre sponding  to  solutions  M 
«nd  [BJ-  In  particular,  represents  a force  doublet  along  the 

s-axis  together  with  a center  of  compression;  both  located  at  the  origin, 
le  note  that  solutions  M (n—  0)  lj  2,  • • • )|  and  pg  (r  = 1,2, 
are  self-equilibrated  in  the  sense  that  the  associated  tractions  on  any 
closed  surface  surrounding  the  origin  are  statically  equivalent  to  null. 

In  the  specific  problem  to  be  treated  presently  we  need  to  construct 
particular  solutions  which  leave  the  plane  i = 0 (p  = 0)  free  from 
surface  tractions.  To  this  end,  we  recall  that 

P2n  + 1(0)  = °»  P2n(0)  = °*  <n  = °»  1»  2»  •••)•  (16) 


Furthermore,  Equations  (16),  by  virtue  of  the  Legendre  recurrence  re- 
lations. 


(2n  + l)p  Pn 


= (n+  1) 


Pn+  1 + 


nP 


n-  1* 


-2 

P 


nP 


n - 1 


inply 


Y (17) 


7 — — 

See  [_?J,  pp.  106-107. 
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P2n+l'°>  = <2n  + 1J  V0) 


P2n+2^ 


p (n\ 

2n  + § P2n'°^ 


P2n<°> 


v ' ?n  5 

22n(L2J 

— 0,  1,  2,  ♦ . . ) . 


Substitution  of  (16),  (18)  into  (12),  (1U)  thus  yields  for  the  nor 


and  shearing  stresses  of  solutions  M and  £b  J,  on  the  plane  p = 0, 


the  following  values: 


For  solution 


ge  = - (2VnO"  ’'re'0! 

r 


for  solution  C*2„+ll 


-re* 

T 


for  solution  ^2n+ 

t2n^>(^-l*2vJW0)>  r 

T 

for  solution  QB2n  + 23> 

,=0,  -r  . - (2n+l)(2n+3)(lm2+  12n+7+2»>; 

% (2n  7 2)7^ 


F,„( 0).  (22) 


In  view  of  Equations  (19)  to  (22),  it  is  clear  that  the  linear  combina- 
tions. 


P>J  = «2n  l>2  J * <2n  + D CB2„+  J 
I? J = a'23+2  J + (Jn  + LB2n  + 2^ 


(n  = 0,  1,  2,  •••), 


L. 


***  V-  ‘-c*  .n.v 


w.p'T  ■-  -^V 
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where 


o^ns  (2n  + 1)Z  - 2 + 2>>, 


(2U) 


satisfy  Equations  (3),  that  is,  clear  the  plane  z s 0 from  tractions. 
From  Equations  (11),  (12),  (13),  (lb),  we  obtain  for  Q)J, 


2Gu 

r 

= - P2n  ♦ (2»  •>  2)(2»  •*  5 - U V)  P2n+  2] 

20^ 

■ - 7^[^n p2„  + <2"-2+  *•  + 1)  p2n+ *] 

°r  * 

(2n'*rElT  n [“»nP2n  + ?2n  P2n  + 2] 

1 

r 

[«^n+  <2"  + l)(2n-2  + UV)]p2n+1 

~ (2n 

+ !K2"+  2)[«5„  P2n+  <tt2n-!ln+  2 

ii 

Un  + 3 

-7^1 

£(2n  + l)(2n  + 2)(1  - 2 V)  Pjn+  ? 

♦ <2n  - 2*>  P2n+  l] 

7tB~ 

(2n+  2)  *2n  P2n  + *2n  + ^ ®2n+ 1 P2n+  2]’ 

and  for  [ 

20u 

r 

- - 7M  [«*n*2p2«+  1+  <2n+  3)(2n+6-U^)P2n4  3] 

2thk 

3 [«2n+  2P2n+  1+  <2n+  2><2n-  1+  MP,B+  3^ 

(25) 


(26) 


(27) 


11 


°r 


% ' -JZ+T 

T 


= (2n  + 2)(2n  + 3)  L p , « p 1 

in+  U ^ 2n+  2P2n+  1 ?2n  + lP2n+  3j 

*2n+2+(2n+2><2n-l*U>'>]P2n+2 
- (2n+  2)(2n+  3)  [«-2n+ 2P2„.,  ! 

+ (a2n+l-2n',l-ll'')  P2n*3] 

<ry='7^[<2n'1<2'')  P2n+  2 

♦ (1  - 2 jJ)(2n  ♦ 2 )(?n  + 3)  P2n+ 3] 

'/r©  ” r2if+  U *20+2  [<-«>*.♦,  * ^2n+  P2n+  l]» 


(28 'j 


where 


pn  = (n  + 2)(n  ♦ 5)  - 2t>. 


(29) 


Solutions  Q)  ]]  snd  £eJ  are  self-*  4 ..^librated  in  the  upper  half- 
cpace,  in  the  sense  that  the  resultant  of  tne  corresponding  tractions 
vanishes  on  a hemisphere  centered  at  r = 0 and  lying  in  the  region 
z — 0.  To  show  this,  we  observe  that  for  any  axl symmetric  stress  field 
this  resultant  is  a single  farce  coincident  with  the  z-sjcis,  the  magni- 
tude of  which  is  given  by 


Z(r)  = - 27 Xr 


[p°-r  - P KJ,  dp. 


r©j 


Do) 


where  Z(r)  is  positive  if  its  sense  is  that  of  the  positive  z-axis» 
Equation  (30),  after  a trivial  computation,  yields  for  the  first  Bous- 
sinesq  solution,  defined  by  Equation  (5), 


Z(r)  = 2rf  r 


30 

Br 


p=0 


(31) 


whereas  for  the  second  Boussinesq  solution,  defined  by  Equation  (6), 
we  reach 


12 


n1 


Z(r)  = hrr(  1 - t>)r£ 


dp. 


(32) 


Equations  (31),  (32),  together  with  (8),  (9),  (23),  (2i*),  confirm  the 
self-equilibrance  in  the  upper  half-space  of  CEJ-  "•  note  that 

self-equilibrance,  in  the  present  sense,  is  not  a necessary  consequence 
of  the  requirement  that  the  tractions  of  a solution  be  zero  on  the  plane 
p = 0,  with  the  exception  of  the  origin.  The  well  known  solution  of 
Boussinesq  corresponding  to  a concentrated  normal  load  on  the  half -space 
conforms  to  this  requirement  but  is,  evidently,  not  self -equilibrated. 

The  physical  significance  of  the  singular  solutions  Q)  ],  0E  ] 
can  be  exhibited  readily.  For  example,  [D^]  corresponds^  to  two 
mutually  perpendicular  tangential  force  doublets,  acting  at  the  origin, 
on  the  plane  p = 0.  The  physical  significance  of  the  other  members  of 
these  two  aggregates  can,  similarly,  be  established  through  suitable 
limit  processes  applied  to  the  solutions  corresponding  to  various  con- 
figurations of  normal  and  tangential  forces  acting  on  the  plane  boundary 
of  the  half -space. 


8 


0See  [9],  p.  191. 
^See  [10]. 


III.  SOLUTION  OF  PROBIZM 


Ve  assume  the  solution  £s]]  to  the  problem  characterised  by  the 
boundary  conditions  (1),  (3),  (U)  in  the  form. 


m - ra + m- 


Solution  [U]  corresponds  to  the  uniform  stress  field. 


(33) 


°x 


xy 


C3U) 


or,  with  reference  *o  spherical  coordinates,  to 


°r 


rr& 


= pp 


j 


where  p,  p are  defined  in  Equations  (10),  while  is  the  solution 
to  the  "residual"  problem,  the  definition  of  which  is  implicit  in  Equa- 
tion (33).  Thus,  is  characterised  by  the  boundary  conditions  (3), 
in  addition  to  the  requirement. 


rr»  , <Ty,  7^e  o as  r-woo, 
and  must  satisfy 

<rr  = - f,  » - PP. 

for  r=l,  O^pM,  0 - y-  27f. 


(36) 


10 


r (37) 


^The  stresses  7T,r  and 
*trv  and  have  been  omitted. 


7*,^  are  sero  by  virtue  of  rotational  sym- 

y& 


13 


fPMBHP**?  VWW-r  f.'iy  ^WIWW^PI  HWIBi 


The  tractions  of  £Rj  on  the  surface  of  the  pit  have  zero  resultant. 
Since,  furthermore,  QT]  must  leave  the  surface  1 - r < oo,  p=0  free 
from  tractions,  we  are  led  to  expand  [R^  in  terms  of  the  two  aggre- 
gates of  solutions  CU  W-  Equations  (23)  to  (29),  which  possess 
both  of  these  properties.  Thus,  we  assume 


ra- 


il P>J  + bn  W 


where  an,  (n  = 0,  1,  2,  ...)  are,  as  yet,  arbitrary. 

In  view  of  the  fact  that  each  (n  = 0,  1,  2,  ...)  satis- 

fies the  boundary  conditions  (3),  (36),  the  same  is  true  of  QT]  pro- 
vided the  series  in  Equation  (38)  converge  suitably.  With  a view  toward 
meeting  the  boundary  conditions  (37),  we  observe  that  they  may  be  written 


°r  = ” f Po  + f P2»  rre  " " 5 P P2> 


for  r=l,  0-p-l,  0-y-2  7T. 


Furthermore,  the  expansions,' 


P2k+  l(p)  = (2k  + 1}  P2k(0) 


4n)  P2n(p)' 


P2k+l(p)  = (2k+  V'  P2k(0) 


,nv  (Un+l)P  (0) 

= — , o < P * l, 

(2k  ■+  1 - 2n)(2k  + 2 + 2n) 

with  the  aid  of  Equations  (28),  lead  to  the  following  Legendre  series 
representation  of  cr,,  7^  of  W on  r = 1,  0-p-l, 


^ P2n<e>- 


»e,  for  example,  QT],  pp.  38,  39i 
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00 

<rr  = (2k+l)(2k+  3)P21c(0)  ' (2k-*-2)a,2k+2  6)£n) 

a.  [(W15 

rie  - (2kt  1K?k+ »vo>  [4n)-  4*  i] 

for  r = 1,  0-p^l,  0 & y & 27f . 

> 

Application  of  the  boundary  conditions  (37)  to  QQ  in  Equation  (38), 
by  aid  of  Equations  (26)  and  (1*1),  yields  a doubly  infinite  system  of 
equations  for  the  coefficients  of  superposition  an,  bR  (n  * 0,1,2, 

These  equations  are  equivalent  to  the  following  sett 


a-i  5 1 


(2n-2k-l)a.2|[+2  4)'n> 


[(2k*  3)pa+1-  (2n*  • ck 


(n  *=  1,  2,  3,  ...), 


an  C ^7^  f ^ ' [<2k+  2>*2n  -1 


“2np2n-2]Qt2k+2 

+ [2nfen-2flt2k  + 2"(2k+  3)<t2n-l?2k+  lj^k+^l 


(n  0,  1,  2,  •••),  J 

2 — — 

Note  that  although  Equations  (1*0)  hold  only  for  the  open  interval 
0 < pj * 1,  Equations  (1*1)  are  valid  for  the  closed  interval  since  Ti- 
of  [fc,  ~1  vanishes  for  d = 0, 
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g(n) 


where  b'“7  is  the  Kronecker  delta  and 
m 


ck  = (2k  + l)(2k  + 3)  P2k(0)  bk,  = 1, 

An=  (2"-D'2n+J)^n[2"Pin-2- 


(1*3) 


The  elimination  of  a^  between  the  two  systems  of  equations  given 
in  Equations  (1*2)  yields  an  infinite  system  of  equations  for  tho  un- 
knowns ck  (k  = 0,  1,  2,  ...).  Since  these  equations  are  at  once  ob- 
tainable from  Equations  (1*2),  they  may  be  omitted  here.  Once  the  ck 
have  been  determined,  the  coefficients  of  superposition  a^,  bQ  follow 
from  Equations  (1*2),  (1*3).  The  values  of  these  coefficients,  together 
with  Equations  (21*)  to  (29),  (33),  (35),  (38),  constitute  the  complete 


solution  of  the  problem. 


IV,  STRESS  EVALUATIONS,  ACCEIZRATION  OF  CONVERGENCE,  NUUERIC1L  RESULTS 


The  infinite  system  of  equations  for  the  c^,  which  is  immediate 
from  Equations  (U2),  was  dealt  with  by  the.  usual  segmentation  process, 
for  ■>>  = Thus,  systems  of  n equations  in  n unknowns,  which  cor- 
respond to  principle  sub-matrices  of  rank  n,  were  solved  for  values  of 
n ranging  from  10  to  26.  The  systems  of  equations  for  n = 25  and 
n = 26  yield  coefficients  c^  which  agree  to  the  number  of  figures 
shown  in  Table  1.  These  values  were  used  in  the  subsequent  stress  com- 
putations. 


Table  1.  Numerical  Values  of  c for  l/U. 

n 


n 

c 

n 

0 

- 0.158,199,1 

1 

O.OU6,7U5,3 

2 

0.006,036,9 

3 

0.001,551,8 

h 

O.OCC,U91,7 

5 

0.000,163,5 

6 

0.000,01:7,3 

The  series,  Equations  (U2),  for  a^,  converge  slowly  and,  conse- 
quently, the  same  applies  to  the  series  for  the  individual  stress  com- 
ponents. For  the  purpose  of  accelerating  the  convergence,  the  following 
method  of  stress  computation  was  adopted:  Let  sn(r,  ©)  be  a repre- 

sentative stress  component  of  W-  If  the  first  of  Equations  (U2)  is 
written  in  the  form, 

oo 

an  _ y “ + ^ T ^ cj{»  2,  3,  •••)>  (Ui) 

\r  E Ct 

••  W 


17 


16 


the  corresponding  representative  stress  component  s(r,  &)  of 


oo 


" an - 1 CDn- J becom8S 

n = i 


oo 


X - n-l  n- X 


n = 1 
oo 


oo  oo  r*  oo  -i 

Vra  1 + y y g<k> . 1 


m 


The  slow  convergence  of  the  double  series  in  Equations  (U5)  is  readily 
traced  to  the  slow  convergence  of  the  series. 


oo 


= gnk>  W 

00 


(U  6) 


Since  the  general  term  in  the  series  for  qv  ' is  available  from  Equa- 
tions (26),  (U2),  (UU),  Hummer' s transformation^  can  be  used  to  acceler- 
ate  the  convergence.  Thus,  a sequence  of  functions  h^  '(r,  Q)  is  de- 


termined such  that 


g(k)  s 

lim  n rrr~~  = 1 , (k  - 0,  1,  2,  ...),  (1*7) 


n — ► oo 


and 


H 


00 


oo 


yvk> 

£1  n 


(U8) 


admits  a closed  representation.  One  then  has 

qw  . „«  + rg«  , . h(«i 

X y l_  n n-l  n J 
n = 1 


(U9) 


^ *> 

‘‘"''See,  for  example,  £ll]]>  P»  2U7 


The  series  in  Equation  (1x9)  converges  faster  than  that  in  Equation  (1*6). 

K 

We  illustrate  the  foregoing  procedure  with  the  stress  component  0 
on  the  axis  of  symmetry,  which  is  of  primary  interest  in  the  present  in- 
vestigation. Thus,  let  s in  the  preceding  discussion  be  <Sq  for 
p = 1.  Since 


Pn(l)  = 1,  p^(l)  = J^-i^,  (n  = 0,  1,  2,  ...), 

we  have  from  Equations  (26),  (1*2),  (1*3),  and  (50),  with  y>  - l/l*, 

(k)  n(Un-3)(l6n2-l)(2n  £k+  V P2n(0) 

8n  n_1  ( l6n2  + i*n+  3)(2n+  2k+  2)(*n- 2k- 3)(2n+  2k+  l*)r2n+  * 

4=<2k+3>P2k+1-(2k+2^2k+2 

■ 

^ = (2k  + 3)(2k  + 2)  p2k  + 1 - 2(2k2+8k+7)ct2k+2 


(n  — 1,  2,  3,  • ••),  (k  0,  1,  2,  ...)• 


for  the  general  term  in  Equations  (1*6).  We  now  choose  h'  ' in  Equa- 

n 

tion  (1*7)  in  the  form, 


(k)  _ CUn2  tk  * 2^  * V P2n(0) 
h"  <2n  + !)(,„ * 

(n  = 1,  2,  3,  • ••),  (k  = 0,  1,  2,  ••«), 


where  1^,  are  independent  of  n.  Equation  (52)  evidently  conforms 

to  the  limit  condition  (1*7)  for  arbitrary  choice  of  the  sequences  1^, 

(k) 

and,  as  will  be  seen  later,  the  corresponding  Hv  7 in  Equations  (1*8) 
admits  a closed  representation.  It  is  apparent  from  Equation  (1*7)  that 


the  best  possible  convergence,  within  tho  present  choice  of  h^,  of  the 
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series  in  Equation  (U6),  is  secured  by  determining  L^,  in  such  a way 
that  the  coefficients  of  the  highest  power  of  n in  the  numerator  of 


(k)  . (k) 

g s , - h'  ' 
6n  n - 1 n 


(53) 


n 


vanish.  This  requirement  leads  to 


4k-  (2k^  2)  tk 

1^  = (16k2  + liUk  + 2U)  ik  - (2k  + 2)  ik 


(lc  0j  ly  2 , •«•)• 


J 


(5U) 


,0c) 


We  now  proceed  to  establish  H ' in  closed  form.  To  this  end  we 


recall1"’  that 


oo 


yr  a2"  p2n<o) 


^ ^ 


Vi+ 


- 1,  0 — JJ. 


* 1. 


(55) 


H- 


Two  successive  integrations  with  respect  to  jx  of  Equation  (55),  and 

2 

subsequent  division  by  /JL  > yield 
- A2"  P2n<°>  — -1 


00  ,,2n 

y-JL 

(2n  + 


sinh  " /uL  _ vr+ 


n = 1 


(2n+  1)  (2n+  2) 


A 

♦Vr 

/? 


(56) 


Similarly  we  obtain,  by  successive  differentiation  with  respect  to  jju 
of  Equation  (56), 


^See  [V],  p.  15. 
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y-  Vvoi 

(2n+  l)(2n  + 2) 
(2n)2^2n  P2n(0) 


Binh~^/x  + 2V1  + 


A 


X A 


05 


2 


f=-*  (2n  + l)(2n  + 2) 


= Binh~1/££ . (** 4 i/sil.  + J» 


In  view  of  Equations  (I48),  (52),  H 


Z4 

00 


M- 


2VI  + 


Z4 


M- 


(57) 


in  the  present  instance  is  ob- 
-1 


tained  by  replacing  ^66  in  Equations  (56),  (57)  with  r , subsequently 
dividing  these  equations  by  r^,  and  forming  the  linear  combination 


apparent  from  Equation  (52)«  This  computation  results  in 

!k 



‘ r2V7T7 


H(k)  = < \ * V 


sinh"1  (i) 


+ (U  - 2I1t  + l'k)(7  - 


(k  — C,  1,  ?,  ...)• 


(58) 


To  demonstrate  the  improvement  in  the  convergence  of  the  original 
(k) 

series  for  q , Equations  (U6),  in  the  present  instance,  we  observe 

2n+  3 3/2 

that  the  general  term  in  this  series  is  of  the  order  (r  n ' ) 

(k)  (k) 

as  n —*■  ao,  whereas  g^  ' 8n - 1 "*  ^n  ‘ in  Equations  (1*9)  here  is  of 

..  - , 2n+3  9/2 x"1 

the  order  (r  n ) . 

An  analogous  procedure  was  used  in  the  determination  of  0"  for 

r 

r = 1,  p — 1,  which  was  confuted  merely  as  a check  on  the  satisfaction 
of  the  boundary  conditions.  Within  the  nurJoer  of  terms  considered  in 
the  numerical  evaluations  of  the  series  solution  established  earlier, 
this  value  of  <r.,  which  should  be  zero,  was  found  to  be  -It  x 10”^  for 
V =^.  This  indicates  satisfaction  of  the  boundary  conditions  to  at 
least  four  significant  figures.  The  variation  of  the  maximum  stress 
along  the  axis  of  symmetry  is  shown  in  Table  2.  These  results  are 
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pared  in  Figure  2 with  the  corresponding  stress  values  given  by  Maun sell 
£l[]  for  the  plane  analog  of  the  present  space  problem.  As  usual,  the 
stress  concentration  in  the  space  problem  is  found  to  be  less  severe  in 
intensity  and  more  localized  in  character  than  it  is  in  its  two-dimensional 
counterpart. 


Tabla  2.  The  Stress 

for  p = 1,  **  = l/U* 

r 

% 

1 

2.23 

1.15 

1.63 

1.35 

1.29 

1.55 

1.15 

2 

1.0U 

3 

1.00 

It  is  interesting  to  note  that  the  maxi  mum  value  2.23  of  ^ at 
the  bottom  of  the  pit  is  7 per  cent  larger  than  the  maximum  stress  con- 
centration induced  by  a spherical  cavity  in  an  infinite  medium  under  the 
same  loading  conditions,1^  and  for  y>  = l/Ii.  In  contrast,  the  maximum 
stress  concentration  at  the  bottom  of  a semi-circular  notch  in  a plate 
under  unit  tension  at  infinity  is  only  2 per  cent  larger  than  the  coi>- 
responding  value  appropriate  to  the  infinite  plate  with  a circular 
hole.17 


16See  [13],  pp.  359-362. 
17See  [1],  [2],  [3]. 
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